The Origin of Magnetic Interactions in Ca3Co206 
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We investigate the microscopic origin of the ferromagnetic and antiferromagnetic spin exchange 
coupUngs in the quasi one-dimensional cobalt compound Ca3Co2 06. In particular, we establish a 
local model which stabilizes a ferromagnetic alignment of the S — 2 spins on the cobalt sites with 
trigonal prismatic symmetry, for a sufficiently strong Hund's rule coupling on the cobalt ions. The 
exchange is mediated through a 5* = cobalt ion at the octahedral sites of the chain structure. 
We present a strong coupling evaluation of the Heisenberg coupling between the 5 = 2 Co spins 
on a separate chain. The chains are coupled antiferromagnetically through super-superexchange via 
short 0-0 bonds. 



Recently there has been renewed interest in systems 
exhibiting magnetization steps. In classical systems such 
as CsCoBrs one single plateau is typically observed in the 
magnetization versus field curve at one third of the mag- 
netization at saturationi This phenomenon attracted 
considerable attention, and Oshikawa, Yamanaka and 
Affleck demonstrated that Heisenberg antiferromagnetic 
chains exhibit such magnetization plateaus when embed- 
ded in a magnetic fieldn^ These steps are expected when 
Nc{S — m) is an integer, where Nc is the number of sites 
in the magnetic unit cell, S the spin quantum number, 
and m the average magnetization per spin, which we shall 
refer to as the OYA criterion. The steps can be stable 
when chains are coupled, for instance in a ladder geome- 
try. In that case the magnetic frustration is an important 
ingredient to their stabilityn^ Plateaus according to the 
OYA criterion are also anticipated for general configu- 
rations, provided gapless excitations do not destabilize 
them^ Indeed several systems exhibiting magnetization 
steps are now knownj^ they all obey the OYA criterion, 
they are usually far from exhausting all the possible m 
values, they all are frustrated systems, and they all can 
be described by an antiferromagnetic Heisenberg model. 
Related behavior has been recently found in other sys- 
tems. For example, up to five plateaus in the magnetiza- 
tion vs. field curve have been observed in Ca3Co206 at 
low temperatur o^'^'^ . However there is to date no micro- 
scopic explanation to this phenomenon, even though the 
location of the plateaus is in agreement with the OYA 
criterion. 

Ca3Co206 belongs to the wide family of compounds 
A'aABOe, and its structure belongs to the space group 
R3c. It consists of infinite chains formed by alternat- 
ing face sharing AOg trigonal prisms and BOg octahedra 
— where Co atoms occupy both A and B sites. Each 
chain is surrounded by six chains separated by Ca atoms. 
As a result a Co ion has two neighboring Co ions on 
the same chain, at a distance of 2.59 A, and twelve Co 
neighbors on the neighboring chains at distances 7.53 A 
(cf. Fig.lSJ'^ Concerning the magnetic structure, the ex- 
periment points toward a ferromagnetic ordering of the 
magnetic Co ions along the chains, together with antifer- 



romagnetic correlations in the buckling a-b planed The 
transition into the ordered state is reflected by a cusp-like 
singularity in the specific heat at 25 Kjii — at the tem- 
perature where a strong increase of the magnetic suscep- 
tibility is observed. Here we note that it is particularly 
intriguing to find magnetization steps in a system where 
the dominant interaction is ferromagnetic. 

In order to determine the effective magnetic Hamil- 
tonian of a particular compound one typically uses 
the Kanamori-Goodenough-Anderson (KGA) rulesi^. 
Knowledge of the ionic configuration of each ion allows to 
estimate the various magnetic couplings. When applying 
this program to Ca3Co206 one faces a series of difficul- 
ties specifically when one tries to reconcile the neutron 
scattering measurements that each second Co ion is non- 
magnetic. Even the assumption that every other Co ion 
is in a high spin state does not settle the intricacies re- 
lated to the magnetic properties; one still has to challenge 
issues such as: i) what are the ionization degrees of the 
Co ions? ii) how is an electron transfered from one cobalt 
ion to a second? iii) which of the magnetic Co ions are 
magnetically coupled? iv) which mechanism generates a 
ferromagnetic coupling along the chains? 

These questions are only partially resolved by ab initio 
calculations. In particular, one obtains that both Co 
ions are in 3+ configurationsii^ Moreover both Co-0 and 
direct Co-Co hybridizations are unusually large, and low 
spin and high spin configurations for the Co ions along 
the chains alternate jJ^ 

Our publication addresses the magnetic couplings, and 
in particular the microscopic origin of the ferromagnetic 
coupling of two Co ions through a non-magnetic Co 
ion. In view of the plethoric variety of iso-structural 
compoundsji^ the presented mechanism is expected to 
apply to many of these systems. We now derive the mag- 
netic inter-Co coupling for Ca3Co206 from microscopic 
considerations. The high-spin low-spin scenario confronts 
us with the question of how a ferromagnetic coupling can 
establish itself, taking into account that the high spin Co 
ions are separated by over 5 A, linked via a non-magnetic 
Co and several oxygens. 

Let us first focus on the Co-atoms in a single Co-0 
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chain of Ca3Co206- As mentioned above the surround- 
ing oxygens form two different environments in an alter- 
nating pattern. We denote the Co ion in the center of the 
oxygen-octahedron Col, and the Co ion in the trigonal 
prisms Co2. The variation in the oxygen-environment 
leads to three important effects. First, there is a dif- 
ference in the strength of the crystal field splitting, be- 
ing larger in the octahedral environment. As a result 
Col is in the low spin state and Co2 in the high spin 
state. Second, the local energy levels are in a different 
sequence. For the octahedral environment we find the 
familiar t2g-~eg splitting, provided the axes of the local 
reference frame point towards the surrounding oxygens. 
The trigonal prismatic environment accounts for a dif- 
ferent set of energy levels. For this local symmetry one 
expects a level scheme with d32,2_^2 as lowest level, fol- 
lowed by two twofold degenerate pairs dxy, dx2-y2 and 
dxz, dyz- However, our LDA calculationai4 show that 
the (i3j,2_^2 level is actually slightly above the first pair 
of levels. Having clarified the sequence of the energy lev- 
els, we now turn to the microscopic processes which link 
the Co ions. Two mechanisms may be competing: ei- 
ther the coupling involves the intermediate oxygens, or 
direct Co-Co overlap is more important. Relying on elec- 
tronic structure calculations, we may safely assume that 
the direct Co-Co overlap dominatesii^ The identification 
of the contributing orbitals is more involved. Following 
Slater and Koster— one finds that only the Sz^-r^ orbitals 
along the chains have significant overlap. However, we 
still have to relate the Koster-Slater coefficients and the 
coefficients for the rotated frame since the natural refer- 
ence frames for Col and Co2 differ. On the Co2 atoms 
with the triangular prismatic environment the z-axis is 
clearly defined along the chain direction, and we choose 
the X direction to point toward one oxygen. This defines 
a reference frame S. The x and y directions are arbitrary 
and irrelevant to our considerations. The octahedral en- 
vironment surrounding the Col atoms defines the natural 
coordinate system, which we call S' . By rotating S' onto 
S one obtains the iz^-r^ orbital in the reference frame 
S as an equally weighted sum of x'y' , x' z' , y' z' orbitals 
in S' . The above observation that the only significant 
overlap is due to the Sz^-r^ orbitals on both Co ions now 
translates into an overlap of the 3z^-r^ orbital on high 
spin cobalt with all t2g orbitals on low spin cobalt. 

We proceed with a strong coupling expansion to iden- 
tify the magnetic coupling along the chain. This amounts 
to determine the difference in energy, between the ferro- 
magnetic and antiferromagnetic configurations, to fourth 
order in the hopping, since this is the leading order to the 
magnetic interaction between the high spin Co ions. As 
explained above we only have to take into account the 
3z^-r^ level on Co2 and the t2g levels on Col. In an ionic 
picture all t2g levels on Col are filled while the d^z'i-r'^ 
level on Co2 is half-filled and we therefore consider hop- 
ping processes from the former to the latter. In the fer- 
romagnetic configuration we include processes where two 
down spin electrons hop from Col to both neighboring 
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FIG. 1: Typical hopping paths for a) ferromagnetic and b) an- 
tiferromagnetic ordering. The displayed ferromagnetic path 
is the only one for ferromagnetic ordering and has the high- 
est multiplicity of all, ferromagnetic and antiferromagnetic. 
There are similar paths for antiferromagnetic ordering but 
with a Hund's rule penalty and lower multiplicity. The path 
in (b) is unique for the antiferromagnetic case and has low 
energy but also low multiplicity. 



Co2 and back again as displayed in Fig. There are in 
total 3x2x2x2x2 = 48 such processes. The inter- 
mediate spin state for Col is in agreement with Hund's 
rule. The energy gain per path is given by: 
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where Acoi and Aco2 denote the crystal field splittings 
on Col and Co2, respectively. The Hund's coupling is Jh, 
assumed to be identical on both. Col and Co2, and U de- 
notes the local Coulomb repulsion. There are no further 
paths in this configuration, besides the one which twice 
iterates second order processes. In the antiferromagnetic 
case the situation is slightly more involved. Here three 
different classes of paths have to be distinguished. The 
first class, denoted al in the following, consists of hop- 
ping events of one up spin and one down spin electron 
from the same Col level. (There arc 3x2x2 = 12 
such paths). The second class (a2) consists of hopping 
events of one down spin and one up spin electron from 
different Col levels (There are 3x2x2x2 = 24 such 
paths). The third class (a3), shown in Fig. QJd, consists 
of hopping processes where one electron is hopping from 
Col to Co2 and then another electron is hopping from 
the other Co2 to the same Col and back again (There 
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are 3x2 = 6 such paths). In total this sums up to 
42 paths in the antiferromagnetic configuration. Conse- 
quently, we have more ferromagnetic than antiferromag- 
netic exchange paths. However the energy gain depends 
on the path. For the classes al and a2, the intermediate 
Col state violates the Hund's rule, and we identify an 
energy gain per path given by: 
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(3C/-2JH-F(Acf,JH))"' 



(5) 



Here F is a positive function which is smaller than Jjj. 
The expression 3U ~ 2 Jh — is the lowest eigenvalue of 
{i\Hco\j) where the states i and j are all possible states 
on Col consistent with two of the c?-orbitals filled and 
three empty. For the class a3 one observes that one does 
not need to invoke a Col ion with four electrons as an 
intermediate state, in contrast to all other processes we 
considered so far. We find the energy gain as: 



E, 



(U + 2 Jh) 



(6) 



Altogether we obtain the difference in energy gain be- 
tween the ferromagnetic and the antiferromagnetic con- 
figurations as: 



A8Ef - 24S,i - 12.B,2 - GEas (7) 



The dependence of E^ — E^^ on Jh for different val- 
ues of the local interaction U is shown in Fig. |21 Us- 
ing Jh = 0.6 eVfii U = 5.3 eV^ t = 1.5 eVfii 
Acoi = 2.5 eV and Aco2 = 1.5 eV^^ we obtain an esti- 
mate for the Heisenberg exchange coupling (for the Co2 
spin 5 = 2): 



(8) 



which is in reasonable agreement with the experimental 
transition temperature of 25 K. 

In this context one should realize that a one- 
dimensional chain does not support a true phase tran- 
sition into the magnetic state. However, as the length 
L of the chains is finite, a crossover into the ferromag- 
netic state may be observed when the correlation length 
is approximately Li^ 

To emphasize the importance of the chain geometry 
we now briefly discuss the hypothetical case where the 
z-axis of the octahcdra corresponds to that of the prism. 
In this geometry there is only one orbital on each Co 
ion which contributes to the hopping processes. In this 
situation the process favoring ferromagnetism shown in 
Fig.^ does not exist, in contrast to the process a3 shown 
in Fig.^, and the resulting coupling is therefore antifer- 
romagnetic. 

In the large class of known isostructural compoundai^ 
the non-magnetic ion is not necessarily a Co ion. If the 




FIG. 2: Energy gain E — E of a nearest neighbor Co2 
ferromagnetic alignment with respect to an antiferromagnetic 
orientation as a function of the Hund's coupling Jh, for a 
typical parameter set of Hubbard U. The grey shaded area 
indicates the interval of Jh for which a Co2-high-spin Col- 
low-spin configuration can be stabilized for the considered 
crystal field splittings^ 



non-magnetic ion is in a 3cP (or AcP) configuration, the 
above argument applies, and the coupling is antiferro- 
magnetic. If the configuration is 3c?'*, all the discussed 
electronic processes contribute, however with different 
multiplicities. Moreover, additional paths have to be 
considered for the antiferromagnetic case. They repre- 
sent exchange processes through an empty orbital on the 
non-magnetic ion. As a result, the ferromagnetic sce- 
nario has fewer paths than the antiferromagnetic, and 
the coupling becomes antiferromagnetic. Correspond- 
ingly, a ferromagnetic interaction can only occur when 
all three orbitals on the nonmagnetic ion participate in 
the exchange process. Obviously the situation we con- 
sider differs from the standard 180 degree superexchange 
mechanism in many respects. 

With the investigation of the interchain magnetic in- 
teraction one first notices that each magnetic Co ion has 
twelve neighboring Co ions on different chains. How- 
ever, as displayed in Fig. |31 there is an oxygen bridge to 
only six neighbors, one per chain. Here the coupling J^^ 
results from the super-superexchangc mechanism (with 
exchange via two oxygen sites), and it is antiferromag- 
netic. Since the Co-0 hybridization is unusally large in 
this system, we expect the interchain magnetic coupling 
to be sufficiently strong to account for the observed an- 
tiferromagnetic correlations. 

From our previous considerations we introduce the 
minimal magnetic Hamiltonian: 



H + JtJS. ■S,)-DY, si, (9) 

with 



i = ±2d 
otherwise 
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and 



if .7 - ±(a + d), ±(6 + J), ±(c + d) 
otherwise. 



Here we use the site vectors a = a(— 1/2, %/3/2, c/(12a)), 
6 = a(-l/2,-V3/2,c/(12a)), c = a(l, 0, c/(12a)), and 
d = a + & + c where a = 9.06 A and c = 10.37 A are the 
lattice constants of the hexagonal unit cell. The Hamil- 
tonian, Eq. Q , also includes a phenomenological contri- 
bution DS^ which accounts for the anisotropy observed, 
for example, in the magnetic susceptibility^*^ 



o 



o 



o 



or© 



©To 

9 



C O 

O Q 



V 



2 6 




*7' 



antifcrromagnetic interchain coupling. Indeed, the lattice 
structure suggests that this magnetic system is highly 
frustrated, since the chains are arranged on a triangular 
lattice. However investigating the Hamiltonian Q re- 
veals that the microscopic mechanism leading to frustra- 
tion is more complex. It is visualized when wc consider 
a closed path TLRT'T — where the sites T and T are 
next nearest-neighbor Co2 sites on the same chain and 
the sequence of sites TLR is located on a triangle of near- 
est neighbor chains. One advances from T to L to i? to 
T' on a helical pathSi formed by the oxygen bridges from 
Fig. 13 Since the structure imposes T and T' to be next 
nearest- neighbors, the frustration occurs independently 
of the sign of the intrachain coupling. 

In summary we established the magnetic interactions 
in an effective magnetic Hamiltonian for Ca3Co2 06. It is 
a spin-2 Hamiltonian, with antifcrromagnetic interchain 
coupling, and ferromagnetic intrachain interactions. The 
latter is obtained from the evaluation of all spin exchange 
paths between two high-spin Co2 sites through an inter- 
mediary low-spin Col site. This mechanism is particular 
to the geometry of the system as is the microscopic mech- 
anism which leads to magnetic frustration. We expect 
that the discussed microscopic mechanisms also apply to 
other isostructural compounds, such as CasCoRhOg and 
CagCoIrOe. 



FIG. 3: Projection of Ca3Co206 in the (100) plane, with 
those oxygen atoms which are closest to the plane. Shad- 
ows indicate the shortest O-O bonds, along which super- 
superexchange processes take place. Yellow (red) large 
spheres denote Col (Co2) atoms, dark (light) small blue 
spheres oxygen atoms located above (below) the Co plane. 



The stability of the magnetization steps results from 
the magnetic frustration which is introduced through the 
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The positions of the considered sites are: T is an arbitrary T' — R + a + d. 

Co2 site on any chain, L — T + b + d, R — L + c + d and 



